Experimental Phonon Band Structure of Graphene using C 12 and C 13 Isotopes 
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Using very uniform large scale chemical vapor deposition grown graphene transferred onto silicon, 
we were able to identify 15 distinct Raman lines associated with graphene monolayers. This was 
possible thanks to a combination of different carbon isotopes and different Raman laser energies 
and extensive averaging without increasing the laser power. This allowed us to obtain a detailed 
experimental phonon dispersion relation for many points in the Brillouin zone. We further identified 
a D+D' peak corresponding to a double phonon process involving both an inter- and intra-valley 
phonon. 

PACS numbers: 81.05.ue, 63.22.Rc, 78.30.-j, 81.15.Gh, 31,30.Gs 



Graphene has attracted a considerable amount of at- 
tention due to the ease in isolating a single sheet of 
graphite via mechanical exfoliation [H, Despite the 
fact that it is one atom thick, exfoliated graphene has 
shown extraordinary electronic, vibrational and optical 
properties that can be used as a novel material for many 
potential applications. Its unique electronic band struc- 
ture constitutes its most noteworthy property, with the 
existence of two degenerate Dirac cones Q, which leads 
to two degenerate valleys (K and K'). Similarly, its vi- 
brational properties, characterized by its phonon band 
structure is also composed of two degenerate Dirac cones 
for the out-of-plane modes. These out-of-plane modes, 
however, do not play any significant role in relation to 
graphene's electronic properties, in contrast to the in- 
plane modes, which can couple to the electronic modes. 
The four in-plane modes are composed of the trans- 
verse (TA) and longitudinal (LA) acoustic modes and the 
transverse (TO) and longitudinal (LO) optical modes, 
which are responsible, for example, for the very high ther- 
mal conductivities observed in graphene [H . 

The most striking example of the interplay between op- 
tical, vibrational and electronic properties, can be found 
in inelastic light scattering (Raman). The photo-excited 
electrons, which live in the Dirac cones, can inelastically 
scatter with the vibrational modes and elastically with 
defects, either within one valley or between valleys and 
back. This leads to a thorough probing of the available 
phase space, for both electrons and phonons. The goal in 
this letter is to use this rich photon-phonon-electron in- 
terplay, in order to deduce the detailed phonon properties 
of graphene, including the in-plane phonon band struc- 
ture, using the well characterized photon and electron 
states. While for graphite, other methods exist to exper- 
imentally characterize the phonon band structure, such 
as inelastic X-ray scattering Q, they do not give enough 
signal with graphene, because of the two-dimensional na- 
ture. Hence, only a very limited number of experimental 
data points for the phonon band structure were obtained 
using Raman spectroscopy on exfoliated graphene 

In order to gain a more complete picture with sufficient 



precision and resolution it is important to avoid finite size 
effects with different edge configurations, which can in- 
fluence both the electronic and vibrational properties of 
graphene. Moreover, the substrate also plays an impor- 
tant role in some of the properties. Hence in order to 
both eliminate substrate effects and to probe large areas, 
we undertook to study Raman scattering for large scale 
chemical vapor deposition (CVD) grown graphene using 
two different isotopes (C 12 and C 13 ) so that we can effec- 
tively exclude and subtract the substrate contributions, 
since a heavier mass downshifts only the vibrational prop- 
erties, while keeping all other properties the same. The 
combination of C 12 and C 13 was also used in CVD- grown 
graphene to investigate the growth mechanism on Ni and 
Cu pj and to study doping effects when prepared as a bi- 
layer graphene [8|. 

For this work we grew graphene monolayers by CVD 
of hydrocarbons on 25 /zm-thick commercial Cu foils. 
The Cu foil is first acid-treated for 10 mins using acetic 
acid and then washed thoroughly with de-ionized water. 
Graphene growth is realized in conditions similar to Li 
and others 0-ll|, but using a vertical quartz tube. We 
used two isotopic methane sources to grow our graphene: 
a 99.99% pure C 12 -methane and a 99.9% C 13 -methane 
(CLM-3590-1) from Cambridge Isotopes Laboratories, 
Inc. Each gas was used to grow graphene on a copper 
substrate. The graphene is grown at 1025°C in 0.5 Torr, 
with a 4 seem H2 flow and a 40 seem CH4 flow for 30 min- 
utes. The methane flow is stopped while the hydrogen 
flow is kept on during the cooling process. The graphene 
film was transferred onto n ++ -doped Si wafer with an ox- 
ide thickness of 285nm. This is realized by spin coating 
the graphene-Cu with PMMA, followed by an etching of 
the Cu foil in an oxidizing solution of 0.1M ammonium 
persulfate((NH 4 ) 2 S 2 8 ). 

Here we performed Raman spectroscopy experiments 
using different laser energies. While there are both 
Stokes and anti-Stokes Raman shifts, we only focused on 
the Stokes lines in this letter, since they are more pro- 
nounced. Hence, we will only consider processes where 
phonons are emitted by electron-phonon scattering. 
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FIG. 1: Photon counts averaged over about 60 Raman spectra 
at different locations of the sample as a function of the Raman 
shift for C 13 -graphene and C 12 -graphene on Si02 covered Si 
substrates. The peaks labeled Si are attributed to the Si/Si02 
substrate. 
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FIG. 2: Difference between the average photon counts for C 12 - 
graphene and C 13 -graphene as an effective way to differentiate 
substrate and surface effects from the intrinsic graphene Ra- 
man spectrum. The red line is a smoothed average of the data 
(black line). 

The Raman spectrum for the C 13 -graphene sample is 
expected to shift downward from the C 12 -graphene by 
a factor of \/l3/12 (the square root of the ratio of the 
atomic masses). The Raman spectra were obtained by 
taking the spectra over many different locations (between 
50-70) of the same sample (of cm 2 size). This allows us to 
confirm its uniformity and to perform a large configura- 
tional average to better resolve very weak peaks, such as 
the D'+D 3 , D'+D 4 and D+D 5 Raman lines, which were 



observed recently 12|. In addition, we identified more 
new weak Raman lines, including D", D 3 , D 4 , D 5 , D 6 , 
D+D' and D'+D 7 , some of which were expected theoret- 
ically ■ The attribution of these new peaks along with 
their peak positions at different laser energies enables us 
to map out a detailed experimental phonon band struc- 
ture. 

Since the substrate can affect the Raman spectrum, 
taking the difference between the spectra of the two iso- 
topes allows to more clearly identify peaks which only 



correspond to graphene even in the presence of substrate 
induced peaks. For the rest of the analysis, which relies 
heavily on the determination of the peak positions, we 
used a combination of methods. For the large peaks, 
such as the 2D, 2D', G and D peaks, we fitted each 
peak after subtracting the substrate spectrum with a 
Lorentzian for every spot on the sample, before deter- 
mining the mean and standard deviation of each peak. 
For very weak peaks this is difficult to do and we then 
only fitted the peaks after averaging over all spots. We 
then compared the peak positions obtained when sub- 
tracting only the average substrate spectrum with the 
peak positions obtained by subtracting the two average 
spectra of the two isotopes. Only if both match, we val- 
idate the peak positions. This is particularly important 
for the graphene peak positions below 1200 cm -1 , where 
the substrate induces large Raman peaks as shown in fig- 
ures 1 1121 It is worth noting that the peaks labeled D 3 , 
D 4 and D 5 in fi gure[2]were not used in the determination 
of the phonon bandstructure discussed below, because of 
the weakness of the peaks. Instead, we used the dou- 
ble resonance peaks labeled D'+D 3 , D'+D 4 , and D+D 5 , 
which are more prevalent. 

To understand the different Raman lines, we can di- 
vide them in three main processes: Process 1 (G-line), 
involves a phonon creation at the T point, with zero mo- 
mentum and hence the electron does not move in the Bril- 
louin zone. This corresponds to the G line at ~ 1586cm -1 
and it is the only Raman process involving a T phonon, 
which leads to a single point on the band structure (BS) 
diagram at the T point in figure [3] 

Processes 2 (intra- valley), are the ones involving a 
phonon emission when the excited electron scatters 
within its valley (K or K'). For processes between oppo- 
site sides of the cones, this leads to a phonon or electron 
wavenumber given by q — ttc/Xvf assuming a linear elec- 
tron dispersion e = 7ii\f|A;| close to the K and K' points, 
where c is the speed of light, vf — 10 6 ms _1 the Fermi 
velocity, A the incident Raman laser wavelength, and k 
the electron wavenumber. These processes can involve 
TA phonons (D 3 Raman line), LA phonons (D 4 Raman 
line), TO phonons (D 7 Raman line) or LO phonons (D' 
Raman line) at momentum q away from the T point (la- 
beled <f* in figured] 

Processes 3 (inter-valley) involve a phonon emission 
for an inter-valley (K to K') scattered electron. Assum- 
ing the electron scatters from opposite sides of the Dirac 
cones (outer processes, which are depicted in the K-M 
diagram) , the corresponding wavenumber of the emitted 
phonon is K-K'+q. If the electron scatters between the 
inner sides (inner processes, which are depicted in the 
T-K diagram), the wavenumber is K-K'-g (labeled ql in 
figured]). Both interactions, K-K'±g lead to Raman lines 
q away from the K-point. These processes can involve 
TA phonons (D 5 Raman line), LA phonons (D" and D 6 
Raman lines), and TO or LO phonons (D Raman line). 



FIG. 3: The top two graphs show the determined phonon 
energies at the corresponding wavenumbers superposed onto 
the theoretical BS taken from [131] for the C 13 (left - scaled 
by yj 12/13) and C 12 (right) samples. In the top left we have 
labeled the 4 in-plane phonon branches, whereas the corre- 
sponding Raman lines are identified in the top right graph. 
The different phonon processes are illustrated in the bottom 
left, where <f* corresponds to an intra-valley phonon and gi to 
an inter-valley phonon within the Brillouin zone depicted in 
the bottom right. 

The D Raman line is usually attributed to the TO branch 

so. 

The various Raman lines can then be mapped onto 
the BS diagram as shown in figure [3J However, some 
of the lines described above are very weak and barely 
visible, even when using the C 13 -C 12 averaging method. 
This is the case for the lines D 3 , D 4 , D 5 , D", and D 7 , 
where we used their double resonance partners instead, 
in order to determine their precise energy. Double reso- 
nance processes, were described by Thomsen and Reich 
for graphite and involve two phonons emitted by the pho- 



ton excited electron 14j. The most dramatic example is 
the well studied 2D line, which is often used to deter- 
mine the number of graphene layers in a given sample 
[i~o| . The corresponding first order D-line is much weaker, 
since it requires the electrons to scatter back across to 
the other valley, which is only possible in the presences 
of short ranged impurities. Hence, while the D-line am- 
plitude increases strongly in the presence of scattering 
impurities, the 2D peak remains the strongest Raman 
line for graphene monolayers. The 2D'-line is very simi- 
lar (but weaker), except that it involves only intra-valley 
scattering. Both the 2D and 2D' processes involve two 
phonons with the same energy but opposite momentum. 
However, there are also a number of phonon pair pro- 
cesses which involve opposite momentum but with dif- 
ferent energies. These are the combination process such 
as D'+D 3 , D'+D 4 , D+D 5 , D+D", and D'+D 7 . They are 
typically stronger than their single phonon counterparts 
because they do not involve an impurity scattering of the 
electron. 



FIG. 4: a) shows the laser excitation dependence of the Ra- 
man lines and their extrapolation to zero energy correspond- 
ing to the high symmetry points of the Brillouin zone. The 
values indicate the energy dependence of the Raman lines in 
units of cm _1 /eV (the errors are in brackets), b) shows the 
equivalence between the energies of the double resonance pro- 
cesses (open squares) and the ones obtained by summing the 
corresponding single phonon processes (full red symbols). The 
missing data points, particularly for 1.9eV, reflects that some 
peaks are too weak for the peak position to be extracted, c) 
illustrates the mechanism attributed to the D+D' line. The 
green crosses on the electronic equi-energy lines (in blue and 
for laser excitations of 1.9, 2.4 and 2.5eV) correspond to the 
possible intra-valley processes, whereas red dots show the pos- 
sible inter- valley processes. 



To show that we can indeed use these double reso- 
nance process in order to determine more precisely the 
Raman lines, we compared the frequencies corresponding 
to (jj£)v + £)n (labeled D"+D'' for a double phonon process) 
with ud" + wop, where D" and are determined by 
the corresponding single phonon processes. We could do 
this for most of the observed double resonance Raman 
lines, which yields ljd»+dp — ^D" + ^d^, within exper- 
imentally accuracy for most lines as shown in figure 0] 
and table |U This justifies the use of the stronger second 
order processes to determine a more precise location of 
the Raman line and map it onto the corresponding single 
phonon BS point in figure [3J A notable exception is the 
line labeled D+D', discussed below, which does not cor- 
respond to processes described above, since it involves 
the combination of an intra-valley phonon (D') and an 
inter- valley phonon (D). 

It turns out that cjd + u>d> > ^Jd+D' > ^>d + ^>g a s 
shown in tablcfl] We can shine light on the D+D' process 
by considering the situation illustrated in figure [4] The 
D'-line (intra-valley) was shown to predominantly occur 
at opposite sides of the Dirac cone along the three K- 
M directions (indicated by small crosses in figure 0] and 
labeled (Jd')- The D-line (inter- valley) was shown to oc- 
cur mainly between opposite ends of the K and K' cones 
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negligible, which is consistent with the theoretically ob- 
tained orientation dependence of only 1 cm -1 13]. This 



TABLE I: List of the difference of some Raman lines for differ- 
ent laser wavelengths and for both isotopes. The errors corre- 
spond to the standard error of the mean except for the lines 
involving the D+D' process at a laser wavelength of 633nm, 
where we could only fit the average over all spots. All num- 
bers are in units of cm -1 . 



along elongated constant energy surfaces [HEH (indi- 
cated by red elongated dots in figure^]) . The wavenumber 
qo connects two of these ends illustrating the outer pro- 
cess. This electron can then either scatter back with the 
help of an impurity to its original momentum (a stan- 
dard D-line process) or it can first emit a r/qo' phonon 
to scatter the electron across the valley before it recom- 
bines through an impurity scattering like the D-line pro- 
cess. The possible values for r] are 0, \ or 1 for scattering 
along the K-M direction (depicted in figure 0] is the case 
1] = h). On average we have r\ ~ i, hence the electron 
loses an energy + uj qD , 

therefore yields the energy 



which we labeled D+D' and 
tabulated in 



OJD, 

table U and corresponding precisely to a; . A simi- 
lar interaction can happen for the inner process, which 
would lead to the same energy loss. The D+D' line is 
therefore consistent with the double phonon process in- 
volving a both an inter- valley phonon and a intra- valley 
phonon. 

In figure |4] we show the dependence of the various dou- 
ble resonance lines on the energy of the laser. This fur- 
ther constrains the assignment of the Raman peaks to 
the different phonon modes. Indeed, The 2D' and D'+D 7 
lines are expected to merge to twice the G-line at the F 
point at zero laser energy, whereas the D'+D 3 , D'+D 4 , 
for example should merge to the G-line at zero energy at 
the r point. This is indeed what happens, further con- 
firming the correct assignments of all the observed peaks. 

Thanks to the high precision in determining peak po- 
sitions of the strong D and 2D peaks, we can see a 
very small difference (but statistically significant) be- 
tween and ujd (~ 2 cm -1 as shown in tablefl] This is 
possibly due to the relative importance of the inner versus 
outer processes, which differ by about 5 cm -1 according 
to Venezuela et al. 13(, since the momentum depen- 
dence of the scattering rates of the electron-phonon pro- 
cess (important for the 2D line) arc different from rates 
of the electron-impurity process only relevant to the D 
line. In contrast, the difference between " 2 „ g/ and lu' d is 



also explains the narrower 2D' peak in comparison to the 
2D peak (see fig. [J). 

Summarizing, we identified several new Raman peaks 
in graphene, using a combination of isotopes and Ra- 
man laser energies on large scale CVD grown graphene, 
which allowed us to do extensive configurational averag- 
ing, without increasing the laser power. In particular, 
we related the peak labeled D+D' to a double phonon 
process involving both an inter- and intra-valley phonon. 
The large number of identified peaks allowed us to map 
out a general in-plane phonon band structure of graphene 
using the excitation energy dependence of the various Ra- 
man lines. 

* Also at the Institute for Microstructural Sciences 
M23A, National Research Council of Canada, Ottawa, 
Ontario K1A 0R6, Canada. 
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